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parameter, namely the perturbation strength could be determined independently from 
level dynamics of the system, thus providing a parameter free agreement between 
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1. Introduction 

The stability of quantum systems under perturbation of the dynamics has attracted a 
considerable amount of attention in recent years. From a practical point of view this was 
motivated by a surge of activity in quantum information, where this question is crucial; 
from a more abstract point of view the implications for "quantum chaos" have been 
studied [1]. Fidelity, as a measure for this stability, has become a standard benchmark 
for the reliability of quantum information processing [2]. Following reference [3], we 
shall consider two unitary operators U'(t) and U(t) for the perturbed and unperturbed 
time evolution, respectively. The fidelity amplitude for some initial state ip(0) is defined 
as 

f(t) = (m\uHt)u'(t)mo)) a) 

and the fidelity by F(t) = \ f(t)\ 2 . In other words, the fidelity amplitude is the cross- 
correlation function of an initial state ip(0) evolving under two different dynamics, 
namely U(t) and U'(t). It can be reinterpreted as the autocorrelation function of the 
same initial state evolving under the echo operator M(t) = U^(t) U'(t) which propagates 
the system forward and backward in time. 

The reason why the fidelity or fidelity amplitude are so attractive is that it gives 
a measure of the average sensitivity to perturbations. In quantum optics there are 
proposals for experiments to measure directly the fidelity amplitude. The idea is to 
study the centre of mass motion of a two level atom in a potential being slightly different 
for the first and the second level. Preparing the atom in a superposition of the two levels 
and measuring the decay of the off-diagonal matrix element of the 2x2 density matrix, 
one obtains the fidelity amplitude of the center of mass motion of the atom in the 
external potential [4, 5, 6]. 

The motivation for the present work comes from recent fidelity studies using 
microwaves [7] and elastic waves [8]. In both types of experiments the waves are coupled 
to the systems by means of antennas or transducers, and reflection or transmission 
measurements are performed. The antennas correspond to external channels, i. e. 
scattering theory is needed to interpret the experiments. It thus becomes evident that 
a generalized concept of fidelity is called for based on the elements of the scattering 
matrix. 

To study the sensitivity of the S-matrix to perturbations we therefore propose to 
define the scattering fidelity of an S-matrix element as the cross-correlation function of 
a matrix element of a perturbed and an unperturbed S-matrix taken in the time-domain 
and normalized by the corresponding autocorrelation functions. We shall see that this 
adequately describes the experimental situation, by presenting an experiment with two 
flat microwave cavities, both chaotic but one with and one without bouncing-ball states. 

We show that scattering fidelity is equivalent to the fidelity amplitude (1) for 
chaotic dynamics and weak coupling; in that case furthermore all scattering fidelities 
are equivalent and we therefore can improve averages by averaging over channels. 
The possibility to replace or improve ensemble averages for S-matrix elements was 
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theoretically discussed in reference [9] and applied in reference [10] to evaluate the 
distribution of absolute values of diagonal and off-diagonal S-matrix elements in a 
problem with many open channels. We then compare the results of our experiment 
with the result for a random matrix model developed in [11]. Indeed the experimental 
results are reproduced without a fit parameter if the somewhat tedious task is carried 
out to determine the strength of the perturbation independently from level dynamics. 

In the next section we define scattering fidelity and its weak coupling limit, then 
we describe the experimental setup (section 3). In section 4 we present a semi-classical 
evaluation of the perturbation strength for our setup, followed by a presentation of 
the experimental results and a comparison with random matrix theory (section 5). 
Particular attention is payed to the effect of bouncing ball states. In the conclusions 
(section 6), we discuss situations where we expect the scattering fidelity to differ from 
the fidelity amplitude, and we consider possible experiments. 



2. Scattering fidelity 

Here, we analyse the stability of a scattering system under perturbation of the 
Hamiltonian. This leads us to the definition of the scattering fidelity for individual 
matrix elements of the scattering matrix. Once the new quantity is defined, we 
ask, under which circumstances that quantity approximates the standard fidelity, 
equation (1). 

Consider a scattering system which can be perturbed in a controlled way, e. g. by 
changing an external parameter. Suppose S'(E) and S(E) are the scattering matrices 
corresponding to the perturbed and unperturbed system, respectively. We denote a 
cross-correlation function of two scattering matrix elements by 

C[Sl b , S'J(E) = (S: b (E>) S' cd (E' + E)) - (S* ab (E')) (S' d (E> + E)) , (2) 

where the brackets denote an energy window and/or an ensemble average. Even though 
the scattering matrix is defined in terms of stationary solutions to a dynamical system, 
we may use its Fourier transform S(t), to obtain a time dependent quantity. Due to 
the convolution theorem [12] such correlation functions are natural candidates for the 
construction of an open systems' analogue to the fidelity amplitude. 

C[S* ab , S'J(t) = JdE e™ Et C[S: b , S'J(E) oc (S ab (ty S' ab (t)) . (3) 

This quantity is dominated by the decay of the autocorrelations, and we therefore include 
a normalization into the definition of the scattering fidelity amplitude: 

e U \ (Sab{t)* S' ab (t)) 

Jab{t) = , • (4) 

V<l-M*)l 2 > (l^(*)l 2 > 

Note that a spectral average taken in equation (2), is equivalent to an appropriate 
smoothing of the time signals (see reference [8]). 

In what follows, we assume that the scattering system can be described by the 
effective Hamiltonian approach [18]. Then, for chaotic systems it holds that the 



Fidelity amplitude of the scattering matrix 



4 



definition (4) approaches the usual fidelity amplitude in the weak-coupling limit. To 
show that, we consider the following parametrization of the scattering matrix 

S ab (E) = 5 ab -i - l — H cS = H int - (i/2) [V V* + T w ] , (5) 

where is the column vector of V corresponding to the scattering channel a. For 
later use, we introduce the absorption width Tw, which is simply a scalar in the present 
context. This is equivalent to model absorption with infinitely many perturbatively 
coupled channels, whose partial widths add up to T w [19]. 

Consider two slightly different Hamiltonians H int = H^Xi) and H' int = if int (A 2 ), 
and the corresponding S-matrices as defined in (5) and denoted by Sat, and S' ab . The 
Fourier transform of S' ab reads: 

S' ab (t) = -i fdE e~ 2niEt \/ (a)t - l — V {b) = -2tt 9(t) V {a)] e' 2 ^'^ V {b) , (6) 

J E ~ H eS 

and similarly for S ab . Here, 9(t) is the Heaviside function. The effective 
Hamiltonians H' eS and if e ff differ only in their Hermitian part which is H' mt or 
.Hint) respectively. Without restricting generality (otherwise, we must perform an 
Engelbrecht-Weidenmuller transformation [13]), we may assume that V = ^2 c y/Wc(v c \, 
where \vi), . . . \vm) are ortho-normal vectors in the Hilbert space of the closed system. 
For the product of two Fourier transformed S-matrix elements, as needed in equation (3), 
we get: 

S ab (tyS' ab (t)=47T 2 e(t)w a w b e- 2 ^ t (v b \e 2 - iH ^ vvU \v a ) (v a \ e -^H Lt t-,vvn ^ _ (?) 

The most unpleasant difference to the usual definition of a fidelity amplitude is the 
projector \v a ) (v a \ between forward and backward propagation. However, under chaotic 
dynamics, and if \v a ) and \v b ) are statistically independent, the energy /ensemble average 
converts that projector into the unit matrix (multiplied by 1/L). We then obtain: 

S ab (ty S' ab (t)) = ^ 9(t) w a w b e~ 2 ^ (v b \ Ul s (t) U' eS {t) \v b ) , (8) 

where U e s(t) and U' cS {t) are the sub-unitary propagators for the respective internal 
Hamiltonians with channel coupling. This looks formally like a fidelity amplitude for 
the effective Hamiltonians of the system. However, even without any perturbation, 
this quantity is not constant, but yields the autocorrelation function (i. e. the power 
spectrum) . This is the motivation to define the scattering fidelity f ab for both diagonal 
and off-diagonal S-matrix elements via the relation 

s ab (ty s' ab (t)) = f ab (t) (\s ab (t)\ 2 ) , (9) 

or, more symmetrically, by equation (4). 

In the limit of weak coupling to the antennas, the scattering fidelity f ab (t) 
approaches the fidelity of the closed system. More rigorously, it must be required that 
the time scale for decay associated to the coupling V is small as compared to the time 
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scale for the decay of the scattering fidelity. In this case the wall absorption, described 
by F w , is the only remaining effect leading to: 

A 2 

~ \ 0(t) w a w b e- 2 ^ /(f) ~ f(t) (\S ab (t)\ 2 ) , (10) 

where f(t) is the fidelity amplitude of the closed system with respect to the initial state 
\vb). This argumentation assumes that \v a ) and \v b ) are statistically independent, in 
particular a ^ b. The result holds, however, for the case a = b as well, provided that 
\v a ) is effectively a random state in the eigenbasis of i?i nt . 

For chaotic dynamics, we have shown that in the limit of small antenna coupling, 
the scattering fidelity f a b(t) approaches the fidelity of the closed system. In this case, 
there is no restriction on the strength of the perturbation. However, an alternative 
derivation may be based on the rescaled Breit-Wigner approximation [20, 19]. That line 
of argument is presented in Appendix A. It is valid in the perturbative regime, only, but 
permits larger antenna couplings (as long as the rescaled Breit-Wigner approximation 
remains valid). The experimental results presented below have been obtained in either 
one of these regimes. 



3. Experimental setup 

Since a detailed description of the general experimental technique can be found e. g. 
in [17], we concentrate on the aspects relevant in the present context. Reflection and 
transmission measurements have been performed in a flat microwave cavity, with top and 
bottom plate parallel to each other. The cavity is quasi-two-dimensional for frequencies 
v < v max = c/(2h), where h is the height of the billiard. In this regime there is 
a complete equivalence between the stationary wave equation and the corresponding 
stationary Schrodinger equation, where the z component of the electric field corresponds 
to the quantum mechanical wave function, 



(A + E) V(x, y) = , E = 




with Dirichlet boundary conditions. 

To describe the experiments adequately, it is essential to take into account that 
the microwave system is in fact an open system. The microwave field in the cavity is 
constantly fed via one antenna, while the second antenna and, most importantly, the 
non-ideal cavity walls act as sinks. Scattering theory, originally developed in nuclear 
physics [18], can be applied directly to open microwave billiards as well [14], if absorption 
is properly taken into account [19]. With H int being the Hamiltonian of the closed 
billiard, and two antennas projecting into the resonator at the positions fi and r 2 , the 
scattering matrix is precisely of the form given in equation (5). In the present context, 
y(a) correS p 0nc l s to the antenna at position f a . For antenna diameters small compared 
to the wavelength the matrix elements Vj a are proportional to ipj(r a ), the wave functions 
of the unperturbed system at the antenna positions. Absorption is included by adding 
the width Tw (as a scalar) to the imaginary part of H e g. 
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B=475mm B=438mm 

Figure 1. Geometry of the billiards. In the right billiard bouncing balls have been 
avoided by inserting additional elements (see text for dimensions). 

Both antennas consisted of copper wires with a diameter of 1 mm, projecting 3 mm 
into the resonator; their positions in the billiard were fixed for all measurements. An 
Agilent 8720ES vector network analyzer was used to determine the complete S-matrix. 
Measurements were taken in the frequency range from 3 to 17 GHz with a resolution of 
0.5 MHz. 

The geometries of the billiards studied in this paper are shown in figure 1. One 
billiard consists of a rectangular cavity of length L = 475 mm, width B = 200 mm and 
height h = 8 mm, a quarter-circle insert of radius R\ = 70 mm, and a half-circle insert 
of radius R2 = 60 mm placed at the lower boundary. The position of the latter was 
changed in steps of 20 mm to get 15 different systems for the ensemble average. The 
perturbation of the system was achieved by moving the left wall in steps of 0.2 mm. For 
each of the 15 positions of the half-circle we performed 11 measurements for different 
perturbation strengths. 

The other billiard consists of the same rectangular cavity, but with length L = 
438 mm. It also shares the quarter-circle insert of radius R\ = 70 mm, and on the lower 
boundary the half-circle insert of radius R2 = 60 mm, which again was moved to realize 
an ensemble of 15 systems. Additional elements were inserted into the billiard to avoid 
bouncing-ball resonances: two half-circle inserts with radius R3 = 30 mm, and a wedge 
on the upper boundary. Again the perturbation of the system was realized by moving 
the left wall in 10 steps of 0.2 mm. 

The classical dynamics for the geometry of the billiards is dominantly chaotic, and 
since these are time-reversal invariant systems, we are going to compare the experimental 
results with random matrix predictions for the Gaussian orthogonal ensemble. 

4. Perturbation parameter for shifting of a billiard wall 

In billiard systems the parameter variation is not due to a change of the Hamiltonian, 
but of the boundary condition. It was shown in chapter 5 of reference [14] that both 
situations are equivalent for the case that the parameter variation in the billiard is due to 
a shift of a straight wall. In this case the matrix element of the equivalent perturbation 
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is given by 

(fll) „.,| l *M^, (12) 

where I is the shift of the wall (in x-direction) and L is the length of the shifted wall. 
It follows for the perturbation strength: 

a 2 = (muj 2 ) 

= ? [ L d Vl [ L dy 2 /^(ft»>»M<>>»>\ /^(°.»> <*■(«»■ »>\ ( 13 ) 
Jo Jo \ dx dx I \ dx dx I 

This is true for n^m, but again the result can be easily generalized to the case n = m. 
The averages can be calculated by means of Berry's conjecture of the superposition of 
random plane waves [21]. Close to a straight wall with Dirichlet boundary conditions, 
this yields 

{lpn(xi,yi)lpn(x2,y2)) 

= \ [Jo (W(^-^) 2 + (yi-y2) 2 ) - Jo (A;v / (^i + ^) 2 + (yi-y 2 ) 2 )] (14) 

for the spacial correlation function, where A is the billiard area. Inserting this into 
equation (13) we obtain for the case n ^ m: 

A 2 = <[(#iW] 2 } 



Ak 2 l 



2/2 1 dyi J ^ 2 (^F [ja%1 ^ 2 ' 



/ dy 1 / dy 2 — [J' (ky 2 )} 
Jo J- yi V2 

pL POO 1 

/ d Vl / dy- [J (y)} 2 

Jo J-00 y 



A 2 

Ak 2 l 2 L 8 
~ A 2 37 

The approximation works well for large wave numbers k. Since the mean level distance 
is normalized to one, we can insert A = An for the area of the billiard, and finally end 
up with: 

A2 = §; kV (15) 

The variance of the diagonal matrix elements, ([(/Ji) nn ] 2 ), yields a value twice as 
large. Exactly the same expression was obtained by Lebceuf and Sieber in a completely 
different approach [22] using periodic orbit theory and the ergodicity assumption. It 
is interesting to note that two seemingly unrelated assumptions, one on the Gaussian 
distribution of wave function amplitudes, and the other on the ergodicity of long periodic 
orbits, yield identical results. 
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Figure 2. (top): Logarithmic plot of the correlation function C[S*]_, S' n ] for the 
billiard with bouncing balls (left), and for the one without bouncing-balls (right). The 
experimental results for the autocorrelation are shown in black, while the correlation 
of perturbed and unperturbed system are shown in grey. The smooth solid curve 
corresponds to the theoretical autocorrelation function, and the dashed curve to the 
product of autocorrelation function and fidelity amplitude. 

(bottom): Logarithmic plot of the corresponding fidelity amplitudes. The smooth 
curve shows the linear-response result. For the billiard without bouncing balls the 
perturbation parameter A was obtained from the variance of the level velocities; in the 
other case it was fitted to the experimental curve. 

5. Experimental results 

5.1. Correlation function and fidelity amplitude 

We start with a discussion of the Fourier transform of the correlation function 
C'[<S'ii, S' n ](t) as given in equation (3). Figure 2 shows a logarithmic plot of C'[<S'* 1 , S' n ](t) 
in comparison to the autocorrelation function ClS^, Su](t). The change of area and 
surface due to the perturbation, i. e. the shift of the billiard wall, was taken into account 
by unfolding the spectra to a mean level distance of one. The frequency window of the 
Fourier transforms was 1 GHz wide, and a Welch filter was applied. The correlation 
functions were averaged over an ensemble of 15 billiard geometries. 

The autocorrelation follows the corresponding theoretical curve nicely, where the 
parameters for the wall absorption and the coupling of the antennas had been determined 
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as described in reference [19]. Minor deviations from the theoretical curve can be 
attributed to the small number of levels in this frequency interval. Furthermore, small 
slits between the inserts of the billiards may act as additional decay channels. This is 
in contrast to the more or less homogenous absorption by the cavity wall, which can be 
described by infinitely many weak channels [19]. 

With increasing time, the correlation function C[S* 1: S' u ](t) deviates more and 
more from the autocorrelation function. This behaviour can be described by a product 
of the autocorrelation function and the fidelity amplitude f(t) of the closed system, thus 
confirming expression (10). 

As an expression for the fidelity amplitude, we use the one derived by Gorin et al. 
[11] based on the linear-response approximation, 



where 6 2 (£) is the two-point form factor. Equation (16) describes the fidelity amplitude 
both in the pertubative regime, where the Gaussian decay is dominant, and in the Fermi 
golden rule regime, which shows an initial exponential decay. 

In the following we shall study the fidelity amplitude directly as the deviation 
from the autocorrelation function as defined in equation (4) to achieve a more direct 
comparison with theory. There are two reasons for dividing by the experimental 
autocorrelation function: there is no need to fit the autocorrelation function, and the 
influence of non-generic features in the Fourier transform is reduced. 

The plots in the lower row of figure 2 show that this procedure works very well. 
For the billiard without bouncing balls and the frequency range shown in this figure, 
the perturbation strength was determined directly from the measured spectra via the 
variance of the level velocities. Thus we can describe the experimental results for the 
fidelity amplitude by the linear-response expression without any free parameter. 

As expected, the billiard with bouncing balls shows systematic deviations from 
the random matrix prediction. Only for small times t we find good agreement. We 
therefore concentrate on the results of the billiard without bouncing balls in the following 
subsections. 

5.2. Agreement with the linear-response prediction 

In our experiment the values for the perturbation parameter A vary from A = 0.01 for 
n — 1 and v = 3. .4 GHz up to A = 0.5 for n — 10 and v — 17 ... 18 GHz. Figure 3 shows 
the fidelity amplitude for four different frequency windows. The perturbation parameter 
A has been fitted to the experimental curves. To improve statistics, experimental results 
for /n, / 22 and fyi have been superimposed. The individual curves for these quantities 
were not discernible within the limit of error. 



/(t) = e- 4A2c W. 



(16) 



For the Gaussian orthogonal ensemble C (t) is given by 




(17) 
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Figure 3. Fidelity amplitude of the billiard without bouncing balls, obtained by 
superimposing results for /n, / 2 2 and /i 2 . The smooth solid line shows the linear- 
response result (16), while the dashed line shows the exact result [15]. The perturbation 
parameter A exp has been fitted to each experimental curve. 



For small perturbation strengths, the linear term in the exponential is still close 
to one and thus we observe essentially a Gaussian decay of the fidelity amplitude, as 
seen in figure 3 for A exp = 0.01. With increasing perturbation strength the linear term 
is getting more pronounced, leading to an exponential decay for small times. However, 
for larger times the Gaussian decay again becomes dominant. 

In the range accessible to our experiment (limited by the small ensemble) we find 
good agreement with the linear-response prediction of the random matrix model. For 
the strongest perturbation strength realized in our experiment, deviations of the linear- 
response result from the exact result for the Gaussian ensemble [15] become noticeable 
only at times where fit) has decayed already below 10 _1 . These deviations could not 
be detected by our experiment. 

In [8] Lobkis and Weaver report on a very similar experiment in elastodynamic 
billiards. There, the principal quantity of interest is called distortion, denoted by D(t). 
In our terminology, it is the logarithm of the scattering fidelity: D(t) = —\nf aa (t). 
Interestingly, the experiment itself is performed in the time domain. A short acoustic 
pulse is coupled into the sample via a pointlike pin transducer, and with the same 
device the time signal corresponding to the response of the system is measured. Cross- 
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correlating two such signals, measured at slightly different temperature allows to extract 
the scattering fidelity, where the perturbation is due to the temperature change. In [16] 
the experimental data from [8] are reanalyzed, showing again an excellent agreement 
with the random matrix prediction, derived in [11]. 

5.3. Scaling behavior of the perturbation strength 

The experimental fidelity amplitude was studied for 10 different shifts of the billiard 
wall (Al = n ■ 0.2mm, with n = 1 . . . 10), and a frequency window of width 1 GHz 
was moved through the spectrum. The perturbation strength A 2 entering the fidelity 
amplitude (16) was fitted to the experimental results. 

Figure 4(a) shows the experimental perturbation strength A 2 xp in dependence of 
the number of steps n for three different frequency regimes. We observe an excellent 
agreement with the scaling A 2 oc n 2 as predicted from equation (15). The experimental 
results for A 2 xp in dependence of the frequency range are shown in figure 4(b). They do 
not look quite as nice as the previous results, but they still confirm the scaling A 2 oc zA 

In view of the good agreement with the predicted scaling of the perturbation 
strength, we can now test the time dependence implied by C(t) with better statistics 
by averaging over all data after rescaling: In figure 5, these averages are plotted against 
47r 2 A 2 C(t) and we see that the linear response solution of the random matrix model [11] 
describes all data very well. However, the slope of the resulting curve is not in accordance 
with the prediction and yields a A 2 xp deviating from the theoretical expectation of 
equation (15) by a factor of 0.36. 

This deviation is caused by the fact that we are far from the semi-classical limit, 
for which the expression for A 2 was derived. This is in accordance with numerical 
calculations for the Sinai billiard done by H. Schanz in Gottingen [23]. In cases where 
we determined the variance of level velocities directly from the measured spectra, we 
found the same deviation from equation (15). This shows that the experiment can be 
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12 3 4 

4-n 2 \ 2 C(t) 

Figure 5. Average of the experimental fidelity amplitude on a rescaled axis x — 
4:TT 2 X 2 C(t). The solid line corresponds to g(x) = exp (—ax) with a = Aj? xp /A 2 = 0.36. 

described in a self-consistent way. 
5-4- Influence of bouncing-ball modes 

To study the frequency dependence of the perturbation strength in some more detail, 
we took a smaller frequency window of 0.5 GHz for the Fourier transform and moved it 
in finer steps through the whole frequency range. 

In Figure 6 the ratio A 2 xp /A 2 is plotted both for the billiard with and without 
bouncing balls, where A 2 is the theoretical value of the perturbation strength according 
to equation (15). 

In the right hand plot, the frequencies of the vertical bouncing ball resonances are 
plotted as vertical lines, revealing the origin of the peaks in the frequency dependence of 
A 2 xp . The vertical bouncing ball states correspond to standing waves between the long 
sides of the billiard and remain nearly unaffected by the perturbation. However, the 
unfolding of each spectrum takes the change of area into account, and thus introduces 
a constant drift to these resonances. This leads to a faster decay of the fidelity, thus 
resulting in the strong peaks visible in the figure. 

Note that the result for the billiard without bouncing balls still shows fluctuations. 
The influence of the eigenmodes in the individual frequency window is thus visible 
for this billiard as well. However, there are no systematic peaks like the ones for the 
bouncing-ball modes. 

6. Conclusions 

It was shown in this paper that the scattering fidelity / &(t) is an easily accessible 
quantity which in the weak-coupling limit approaches the ordinary fidelity amplitude 
f(t). It is stressed that it is not necessary to vary the antenna position for this purpose. 
For integrable systems, however, an average over the antenna positions is indispensable 
to obtain f(t). Therefore a determination of the fidelity f(t) from the scattering fidelity 
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Figure 6. A^ xp /A 2 as a function of frequency v. Results for the billiard without (left) 
and with bouncing balls (right). The dashed lines correspond to the eigenfrequencies 
of the associated bouncing-ball states. 



f a b{t) as illustrated here, is probably not feasable in that case. 

All results of the present paper could be described within the limits of the 
exponentiated linear-response approximation [11], allowing a nice scaling of all data 
onto one single curve. Deviations of the decay constant of the such scaled fidelity from 
the predicted value could be quantitatively traced back to limits of the semiclassical 
approximation. 

A general perturbation of the scattering matrix, is not necessarily restricted to the 
internal dynamics, as it may also change the coupling to the channels. For practical 
reasons, the present experimental setup was not well suited to study such a situation, 
and we are planning a further experiment to study scattering fidelity as a function of 
coupling strength. 
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Appendix A. Scattering fidelity in the rescaled Breit-Wigner approximation 

We assume that the spectrum of the unperturbed closed system (Hamiltonian H int ) 
may be taken from the Gaussian orthogonal ensemble (GOE), and consider averages 
over that ensemble. In the rescaled Breit-Wigner approximation [20], the correlation 
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function (7) reads 

(S ab (ty S' ab (t)) = 4tt 2 0(O (j2(*\V\j) 01 e 2 ^' +i W \j) (j\V\b) 

x (b\V\k') (k'le^'^^lk') (k'\V\a)^ , (A.l) 

where we have used Dirac's notation for the eigenbases £V \j) (j\ and J2 k \k') (k'\ of 
iJint and -£^ nt , respectively. The numbers Tj and 1^ denote the diagonal element of 
the anti-Hermitean part of H cS and H' cS in the respective basis. We will relate this 
correlation function to the fidelity amplitude of the closed system Hi nt , or H' int (when 
perturbed). To this end, we need to restrict ourselves to the perturbative regime. In 
this case, \f) = \j) and 1^- = Tj, so that 

(S ab (ty S> ab (t)) = 47^(0 J2 (K V jb V; b V ka e-< T ^ e -^(K-^ 

jk 

= E (e- 2 ^~ Ek)t ) J2 {K Vjb V k * b V ka e -^ +r ^ e-^m-E^ (A 2) 

k j 

The separation of averages in the second line is possible in the perturbative regime, 
only. There E' k — E k is just the diagonal element of the perturbation matrix and hence 
statistically independent of the unperturbed system. The average (exp[— 1-n\{E' k — E k )t\) 
gives the fidelity amplitude, independent of the index k, while the rest is just the 
autocorrelation function of the S- matrix element S ab . 

(S ab (tys' ab (t)) = f(t) (\S ab (t)\ 2 ) ■ (A.3) 
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